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Abstract 

We investigate the low-temperature properties of the spin-1 antiferromag- 
netic Heisenberg chain with bond-alternation by the quantum Monte Carlo 
method (loop algorithm). The strength of bond-alternation at the gap- 
less point is estimated as S c = 0.2595 ± 0.0005. We confirm numerically 
that the low-temperature properties at the gapless point are consistent with 
field theoretical predictions. The numerical results are compared with those 
of the spin-1/2 antiferromagnetic Heisenberg chain and recent experimental 
results for [{Ni(333-tet)(^-N 3 )}„](C10 4 )„ (333-tet=tetraamine iV,JV'-bis(3- 
aminopropyl)-l,3-propanediamine). 

PACS numbers: 75.40.-s, 75.40. Cx, 75.40.Mg, 75.50.Ee 
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I. INTRODUCTION 



For one- dimensional quantum spin chains, there are some theoretical predictions which 
have been confirmed experimentally. One example is the Haldane conjecture 0. Haldane 
investigated the 0(3) nonlinear a model, and predicted that an excitation gap opens for 
integer spin chains but not for half-odd integer spin chains. This conjecture was confirmed 
by numerical calculations and the Haldane gap was observed experimentally |5|-0. 

Another example of such theoretical predictions is the presence of logarithmic corrections 
at very low temperatures for the spin- 1/2 antiferromagnetic Heisenberg chain (S1/2AH). 
This feature was predicted by the renormalizat ion-group approach and was confirmed 
numerically by Bethe Ansatz |[LG|| . Later, this behavior was observed experimentally fll |. In 
this way, the presence of logarithmic corrections at very low temperatures for the spin-1/2 
antiferromagnetic Heisenberg chain was established. 

There is another interesting prediction which has been observed experimentally quite 
recently fl2" |. This is the so-called Affleck-Haldane conjecture p3] , |T4|| . Affleck extended 
Haldane's argument to the bond-alternating chains, and predicted that there will be IS 
gapless points for the spin-S* antiferromagnetic Heisenberg chain with bond-alternation. 
Later, for the spin-1 case, the central charge c at the gapless point is estimated as 
1 by numerical calculations [|Hi|-TT|. This implies that the low- lying excitations at 



the gapless point are described by the level k = 1 SU(2) Wess-Zumino-Witten (WZW) 
model. As a result, the same qualitative low-temperature properties as those of the spin- 
1/2 antiferromagnetic Heisenberg chain are expected at the gapless point. Recent exper- 
iments for [{Ni(333-tet)(/i-N 3 )} n ](C104) n (333-tet=tetraamine N, A^ / -bis(3-aminopropyl)- 
1,3-propanediamine) show that this compound has a structure which is effectively described 
by the spin-1 antiferromagnetic Heisenberg chain with bond-alternation (S1BA) fl8f , and 
that the behavior of the uniform susceptibility is close to what is expected at the gapless 



point [12fl . Thus, this compound is probably an experimental realization of the gapless 
point of the Affleck-Haldane conjecture. In order to clarify this feature in more detail, it is 
necessary to show explicitly how the uniform susceptibility behaves in the low-temperature 
regime for SIB A at the gapless point. Hence, the purpose of this paper is to clarify the 
range of temperatures for which the field-theoretical prediction is valid, and to what extent 
SIB A can explain the low-temperature properties of this compound. 

For the gapless point of S1BA, Singh and Gelfand applied a series expansion technique, 
and estimated the critical value of the strength of bond-alternation 5 as 5 C = 0.25 ± 0.03 
T9| . Later, Kato and Tanaka applied the density-matrix renormalizat ion group (DMRG) 



method, and obtained clear evidence of the Affleck-Haldane conjecture [IB]. They estimated 
the critical value 5 C more accurately as 5 C = 0.25 ± 0.01. As for the excitation spectrum, 
Yamamoto performed the quantum Monte Carlo simulation (world-line algorithm), and 
obtained the dispersion relation |16||. Totsuka, et.al. investigated the low- lying excitation 



spectrum by exact diagonalization, and analyzed the results using the conformal field theory 
|T7|| . They estimated the critical value of 5 as 5 C = 0.254 ± 0.008. 

Before investigating the low-temperature properties at the gapless point, we have to 
determine the critical value of 5 more accurately. In Sec.|H|, we estimate 5 C from the low- 
temperature behavior of the uniform susceptibility and the staggered susceptibility. In 
Sec.[[V], we analyze the low-temperature properties of SIB A at the gapless point based on 
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the field-theoretical prediction. Comparisons with numerical results of the spin-1/2 antifer- 
romagnetic Heisenberg chain are also made. In Sec.0, the numerical results are compared 
with recent experimental results for [{Ni(333-tet)(/x-N 3 )} rt ](C10 4 ) ri . 

In the present paper, we consider SIB A defined by the following Hamiltonian: 

H = J^[l-(-l) ?; 5]5 i -5 i+1 , (1.1) 

i 

where Si denotes the spin operator at site i with spin one (5 = 1). The length of the chain 
and the temperature are denoted by L and T, respectively. We set J = 1 as the energy unit. 

We use the quantum Monte Carlo (QMC) method (loop algorithm) p0| . The simulations 
have been performed in the grand-canonical ensemble. We have typically run 6 x 10 5 Monte 
Carlo steps for measurements after 6 x 10 4 steps. We have made extrapolation for the Trotter 
slice At (At — > 0) as a + b- At 2 using the data at At ~ 1/6, 1/7 and 1/8. The system-sizes 
we have investigated are L = 32, 64, 96, 192, 320 and 400. We have performed calculations 
up to inverse temperature (3 = 100. 



II. FINITE-TEMPERATURE PROPERTIES OF S1BA 

First, let us look over global features of SIB A for finite temperatures. Figure [l] shows 
the uniform susceptibility x(l = 0;T), the staggered susceptibility x(Q = ft'iT) and the 
staggered structure factor S(q = 7r; T) defined as 

x(? = 0;T) = ^((£S?) 2 > r , (2.1) 

i 

X(q = vr;T) = ^((EM^V^V, (2-2) 

S(q = n;T) = j((j:(-iySl) 2 )T, (2.3) 

where £>£ denotes the ^-component of the spin at site in the (l+l)-dimensional space- 
time, and M is the number of Trotter slices defined as M = (3/ At. Here, (■ • -)t denotes the 
thermal average at temperature T. 

The uniform susceptibility x(Q = 0; T) in the low-temperature regime is fitted well by 
x(q = 0; T) oc ^exp(— A/T), except for 5 = 0.25, as shown in Fig.||. This behavior is 
expected when the dispersion relation is E(q)\ q _>o = aq 2 + A [21], where q is the momentum 



measured from the lowest triplet state and a is a constant. Thus, except for 5 = 0.25, the 
low-energy excitation may be explained by a magnon excitation with a finite excitation gap 



A ||3|j22|. This excitation gap A gradually decreases as 5 approaches 5 C ~ 0.25. As for the 
staggered susceptibility x(q = ^\T) and the staggered structure factor S(q = vr;T), the 
maximum value becomes larger as 5 approaches 5 C ~ 0.25. This behavior is expected from 
the Affleck-Haldane conjecture and is consistent with various numerical results |15|-p!7|,[23| . 
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III. ESTIMATION OF THE CRITICAL POINT 



In this section, we accurately determine the critical value 8 C . One criterion is based 
on the expectation that the uniform susceptibility x(q = 0; T) reaches a maximum at the 
critical point if the temperature is low enough. In Fig.|3](a), we show the (^-dependence of 
the uniform susceptibility. The maximum is located near 8 ~ 0.2595 with little size and 
temperature dependence. 

Another criterion for the critical point 8 C is that the staggered susceptibility x(<? = tt; T) 
is a maximum at the critical point at sufficiently low temperatures |24j[. We show the 8- 
dependence of the staggered susceptibility in Fig.^|(b). Near 8 ~ 0.2595, x(q = ^T) also 
reaches a maximum with little size and temperature dependence. Thus, we estimate the 
critical value of 8 as 8 C = 0.2595 ± 0.0005 [||. 

In order to show the validity of this estimation, we compare the low-temperature behav- 
iors of x(q = 0; T) and x(q = n\ T) at 8 = 0.2595 with those at 8 = 0.250 and 0.255 in Fig.|. 
This figure clearly shows that 8 = 0.2595 is closer to the critical point 8 C than 8 = 0.250 or 
8 = 0.255. 

IV. LOW-TEMPERATURE PROPERTIES AT THE CRITICAL POINT 

In this section, we investigate the low-temperature properties of S1BA at the critical 
point, assuming that it is described by the k — 1 SU(2) WZW model (with a marginally 
irrelevant operator) fllQ"| . 

A. Uniform susceptibility 

Before investigating the low-temperature behavior of the uniform susceptibility of SIB A 
at the gapless point, we briefly review the case of the spin- 1/2 antiferromagnetic Heisenberg 
chain (S1/2AH) as an example which is described by the k = 1 SU(2) WZW model. Figure 
H shows the temperature dependence of the uniform susceptibility x(l — 0; T) for S1/2AH 
|T0| . The dotted line corresponds to the fit assuming eq. (f4~T|) , which is expected from the 
renormalization group approach f^fTO] . 



1 1 r 1 HHTo/T) + 1/2), , ,< , , m9x / s 



The second term in eq. ( |4.1| ) is the leading logarithmic correction term due to the marginally 
irrelevant operator (oc Jl ■ Jr). One of the features due to this logarithmic correction is the 
infinite slope in the low-temperature limit. As a result, naive extrapolation of x(q = 0; T) 
as T — > 0, using the data in the low-temperature regime (0.02 ~ T ~ 0.2), does not coincide 
with the zero-temperature uniform susceptibility x(q = 0;T = 0) = 1/(2ttv) Hl0| . Another 



feature is the existence of an inflection point in the low-temperature regime. For S1/2AH, 



the inflection point is near T = 0.087 [10] 



Let us consider the S1BA case. Figure ^] shows the temperature dependence of the 
uniform susceptibility x(q = 0;T) for S1BA at the critical point (8 = 0.2595). This figure 
supports the existence of the logarithmic correction in the sense that naive extrapolation of 
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x(q = 0; T) as T — > does not coincide with the zero-temperature uniform susceptibility 
x(q = 0; T = 0) = l/(2irv), and that there exists an inflection point near T ~ 0.2. We 
have tried to fit the numerical data as eq.(fO), and estimated T ~ 0.34. This value of T 



is smaller than that of S 1/2 AH by about one order of magnitude |10|. 



B. Staggered susceptibility and staggered structure factor 



Next, we consider the low-temperature behaviors of the staggered susceptibility x(q — 
7r; T) and the staggered structure factor S(q = n; T). For S1/2AH, Starykh et.al. confirmed 
that x(q = tt; T) and S(q = 7r; T) in the low-temperature regime behave as follows [26] : 



X (g = vr;T)ocT- 1 [ln(T x /T)] 1 /2 ! 
S(q = n;T)oc[HT s /T)f 2 . 



(4.2) 
(4.3) 



Here, we consider the S1BA case. Figures [7| and || show the temperature dependence of 
the staggered susceptibility x(q — 7r; T) and the staggered structure factor S(q = it; T) for 
SIB A at the gapless point in the low temperature regime. These figures suggest that the 
low-temperature behaviors of x(q — n ',T) and S(q = vr; T) for S1BA at the critical point 
are qualitatively the same as those of S1/2AH. [For comparison, in the insets of Figs. 0(b) 
and ||(b), we show the S1/2AH case.] Thus, the behaviors of x(q = T) and S(q = 7r; T) 
also support that S1BA at the critical point belongs to the universality class of the k = 1 
SU(2) WZW model. 



V. COMPARISON WITH EXPERIMENTAL DATA 



In this section, numerical results are compared with the experimental data for [{Ni(333- 
tet)(/x-N 3 )} n ](C104)„ Hl2| . Figure | shows the temperature dependence of the uniform sus- 
ceptibility for the powder sample of [{Ni(333-tet)(/x-N 3 )} n ](C104) n . The global feature of 
the experimental data is very similar to that of the numerical results for SIB A at the gapless 
point. This suggests that this compound is effectively described by S1BA near the critical 
point. A little difference in the low temperature regime may be due to small anisotropy 
effects or off-criticality of this compound. 



VI. SUMMARY 



We have reported the numerical results for the spin-1 antiferromagnetic Heisenberg chain 
with bond-alternation (S1BA) obtained by quantum Monte Carlo (loop algorithm). We have 
estimated the strength of bond-alternation at the gapless point as S c = 0.2595 ± 0.0005. At 
the gapless point, we have confirmed that the low-temperature properties are effectively 
described by the k = 1 SU(2) WZW model (with a marginally irrelevant operator). The 
numerical results for S1BA at the critical point well explain recent experimental results for 
[{Ni(333-tet)(/x-N 3 )} n ](C104) n , indicating that this compound is close to the gapless point 
of the Affleck-Haldane conjecture. 
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FIGURES 



FIG. 1. Uniform susceptibility x(q = 0;T) [(a)], staggered susceptibility x(q = k^T) [(b)] and 
the staggered structure factor S(q = vr;T) [(c)] of S1BA as a function of temperature T. The 
symbols are 64-site data. The 96-site date for < 5 < 0.25 (0.25 < 5 < 0.5) are joined by solid 
(dotted) lines. The thick solid line corresponds to 96-site data for 5 = 0.25. 

FIG. 2. Uniform susceptibility x(q = 0; T) of S1BA as a function of inverse temperature (3. 
The solid lines correspond to the fit as x{q = 0; T) oc exp(— A/T). Open and solid symbols 
denote 96-site data. Crosses denote 64-site data. 

FIG. 3. Uniform [(a)] and staggered [(b)] susceptibility of S1BA as a function of 5. The bold 
lines are guides to the eye. 

FIG. 4. Uniform [(a)] and staggered [(b)] susceptibility of S1BA near the gapless point in the 
low-temperature regime. L = 192. 



FIG. 5. Uniform susceptibility x(q = 0;T) of S1/2AH for < T < 0.2. The crosses denote 
the result obtained by Bethe Ansatz cited from ref. The dotted line corresponds to the fit 

assuming eq. (0). We choose T = 2.3. Solid and open symbols denote the data obtained by QMC 
loop algorithm. The inset shows x(q = 0;T") of S1/2AH for < T < 2. The solid line in the inset 
is obtained by Bethe Ansatz fllpf . 



FIG. 6. Uniform susceptibility x(q = 0;T) of S1BA at 5 = 0.2595 for < T < 0.4. The data 



at T = are obtained as x(q = 0; T = 0) = 1/(2ttu), where v = 2.46 ± 0.08 |§ and v = 2.39 g7|. 
The dashed line corresponds to the linear fit using the data for 0.02 < T < 0.2. The inset shows 
X(q = 0;T) for < T < 5. 



FIG. 7. Staggered susceptibility x(q = k^T) of S1BA at 5 = 0.2595 in the low-temperature 
regime, (a) Linear plot and (b) logarithmic plot. The solid line corresponds to the fit assuming 
eq.( |4.2D . We estimate T x ~ 8.9. The inset in (a) shows the inverse of x(q = 7r ;2 n )- The inset in 
(b) shows the logarithmic plot of x(q = for S1/2AH obtained by QMC loop algorithm. We 

estimate T x ~ 9.8 for S1/2AH. 

FIG. 8. Staggered structure factor S(q = vr;T) of S1BA at 6 = 0.2595 in the low-temperature 
regime, (a) Linear plot and (b) logarithmic plot. The solid line corresponds to the fit assuming 
eq.( |4.3D . We estimate Tg ~ 20.6. The inset in (a) shows the inverse of S(q = tt;T). The inset in 
(b) shows the logarithmic plot of S(q = tt;T) for S1/2AH obtained by QMC loop algorithm. We 
estimate T s ~ 21.5 for S1/2AH. 
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FIG. 9. Uniform susceptibility x (q = 0;T) of [{Ni(333-tet)(/i-N 3 )} n ](C10 4 )n (open diamonds) 
]i~2fl and that of S1BA for 5 = 0.2595 (solid diamonds). The data at T = are obtained as 
x (q = 0;T = 0) = 1/(2™), where v = 2.46 ± 0.08 @ and v = 2.39 @. The inset shows 
the blow-up region for the low-temperature regime. We choose the g-v&lue as g = 2.46, and 
J/k B = 86 [K] 
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Fig.|(c) 
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Fig-l(a) 



12 



Fig-l(b) 
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Figi(b) 
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Fig.0(b) 
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